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Abstract 

We investigate solutions of Einstein field equations for the non- 
static spherically symmetric perfect fluid case using different equations 
of state. The properties of an exact spherically symmetric perfect fluid 
solutions are obtained which contain shear. We obtain three different 
solutions out of these one turns out to be an incoherent dust solution 
and the other two are stiff matter solutions. 



1 INTRODUCTION 

There is no shortage of exact solutions of Einstein field equation (EFEs). 
However, because General Relativity is highly non-linear, it is not always 
easy to understand what qualitative features solutions might possess. Differ- 
ent people have been working on the investigation of spherically symmetric 
perfect fluid solutions with shear [1-5]. Nearly all the solutions have been 
obtained by imposing symmetry conditions. It is known that non-linear par- 
tial differential equations admit large classes of solutions, many of which are 
unphysical. 

EFEs for static spherically symmetric distribution of perfect fluid have 
been investigated by many authors using different approaches [6]. One ap- 
proach is to prescribe an equation of state, p = p(p) which relates the energy 
density p and isotropic pressure p. In this paper we shall extend this idea 
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of solving EFEs for non-static spherically symmetric spacetimes. We shall 
examine systematically the field equations for non-static spherically sym- 
metric perfect fluid solutions. We obtain three different solutions which have 
shear non-zero. Since most of the models in the literature [6] for spherically 
symmetric are shear-free. Thus it will be interesting to study the non-static 
solutions which contain shear. 

The breakup of the paper is as follows. In section two we shall write down 
the field equations. In the third section we attempt all the possible solutions 
in three classes using different equations of state. In section four we shall 
evaluate kinematic quantities for the solutions obtained. Finally in the last 
section we shall conclude our discussion. 



2 Field Equations 

The non-static spherically symmetric metric has the form [5] 

ds 2 = e 2v ^dt 2 - e 2X ^dr 2 - R 2 (r, t)dtt 2 , (1) 

where dVl 2 = dQ 2 + sin 2 9 dip 2 . 

For a perfect fluid distribution, the energy-momentum tensor is given by 

T ab = (p + P)u a u b - pg ab , u a u b = l, a, b = 0,1, 2, 3. (2) 

where p and p are the energy density and the pressure of the fluid respectively. 
The four velocity of the fluid has the form u a = (e^ (r '* } , 0, 0, 0). The field 
equations for the metric (1) can be written down [6] 

K pR = e~ 2A ((z/" + v' 2 - u'\')R + R" + R'u - R'X') 

-e- 2u {(\ + \ - \u)R+ R + R\- RO), (5) 

R - Rv- R'\=0, (6) 
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where the dot denotes partial derivative with respect to time l tf and the 
prime indicates partial derivative with respect to the coordinate l r\ The 
spatial coordinate V refers to the comoving radius and 'k' is the gravitational 
constant. The consequences of the energy- momentum conservation = 
are the relations 

p' = -(p + p)v', p=-(p + p)(X+2^). (7) 

We now consider the equation of state 

p=( 7 -l)p, p + p = P1 , 1< 7 <2, (8) 

where 7 is a constant. (The limits on 7 result from the requirement that 
the stresses be pressures rather than tensions and that the speed of sound 
in the fluid be less than the speed of light in vacuum). For 7=1, the 
pressure vanishes, so that the equation of state is that of incoherent dust. 
For 7=3, the equation of state is that of photon gas or a gas of non- 
interacting relativistic particles. For 7 = 2, the equation of state reduces to 
the stiff matter case. 



3 Non-Static Spherically Symmetric Solutions 

To simplify the field equations we solve them for some special cases: 
3.1 THE CASE R=R(t), A = A(t) 

This class is identical with the well-known Kantowski-Sachs class of cosmo- 
logical models [7-9]. We can choose R = t without loss of generality, use the 
equation of state given by Eq.(8) and attempt all the possible solutions. 
When 7=1 

This gives p = and the EFEs. (3-6) will reduce to 

Kp = _L + ? e - (A+ J_), (9) 

= 2t v e - 2u - 1 - e- 2u , (10) 
= -e~ 2v [(\ + X - \v)t+ A - i>\. (11) 
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Eq.(lO) can easily be solved which gives 



v — - In 
2 



t 



c-t 



where c is an arbitrary constant. For this value of v, Eq.(9) gives 

put 3 — c 



2t(c - t) 



Substituting the values of A, A , A, v in Eq.(ll), we have 

3c - At t 2 K 



2t(c-ty 2(t-c) 



Solving this, we get 



^— )5 - sin - ^-)*}^ 3 - t 4 + CiVct 3 - t 4 
c t c 



(12) 



where c\ is an integration constant. By substituting this value of p in A and 
integrating, the value of A will become 



A = In 



C — t 1 . -t i l C — t lCi 

c 2 {l - (— -)*sm \-)* + (^) 5 -t) 

v C v r\j 



, C — t l C\ 



(13) 



where c<i is an integration constant. The energy-momentum conservation 
relations (7) are also satisfied by this solution. 
The resulting spacetime is 



ds 2 



c-t 



-dt 2 - 



c 2 {l- 



,c-t 



,t, 



,C — t. l Ci- 



i 2 



— )*sm x (-)2 + (— — )a— } 

6 C I K 



dr 2 - t 2 dn 2 



(14) 

If we take a special case for which A = constant and 7 = 2(p = p), then 
EFEs (3-6) will give 



v = In 



1 - a 2 t 2 



where a is an arbitrary constant. The energy density and the pressure can 
be evaluated as 

1 



p = p = 



Ka 2 t A 



The corresponding metric will be 



ds 2 



a 2 t 2 
1 - ct 2 t 2 



dt 2 - dr 2 - t 2 dn 2 



(15) 



(16) 
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3.2 THE CASE R=R(t), A = A (r,t) 

This class of solution was examined by Korkina and Martinenko [10,11]. In 
this case EFEs will be the same as Eqs.(3-6) except that now A = A(r, t). We 
solve this system of partial differential equations using Eq.(8). 
When 7 = 2 

In this case p — p, EFEs (3-6) will give 



2 z . 

1 + — f +1, 
z 



V 



p 



1 _ Z cy Z cy Z y - 



where y = e 2v{ - l \ z = e x( - r,t \ From Eqs.(18) and (19), we have 



.J I, 



l 2yt yt 2 



= 0. 



(17) 
(18) 

(19) 
(20) 



This second order non-linear partial differential equation can be solved using 
Herlt's method [5] by choosing 



y(t) 



-2u 



1 



n 2 — 1 

where (3 is an arbitrary constant. Eq.(20) then becomes 

z+Az -Bz = 0, 

where 

r 1 - n(n 2 - l)pt- 2n - 3 



A 



B = 



1 + (3{n 2 - l)t- 2n - 2 ' 
n 2 t- 2 - n(n 2 - l)/?r 2 "- 4 



1 + (3{n 2 - 1)£- 2 ™- 2 
Eq.(22) has the special solution z s = t n . Let us substitute 



(21) 
(22) 

(23) 
(24) 



z = C(r,t) t n 



(25) 
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into Eq.(22), we obtain 

C +[A + 2nr 1 ]C=0, (26) 
This can easily be solved, and the general solution becomes 

z(r,t) = f»{ft(r) / dt }, (27) 

where /3i(r) and ^(r) are arbitrary functions of the variable r. The energy 
density p and pressure p can now be computed easily using Eqs.(17) and 
(18). " 

The corresponding metric will become 

ds 2 = [yit)]' 1 ^ 2 - z 2 (r, t)dr 2 - t 2 dtf, (28) 

where y(t) and z(r,t) are given by Eqs.(21) and (27) respectively. It is to 
be noticed that this solution corresponds to the solution obtained by Herlt 
[5]. 

3.3 THE GENERAL CASE R=R(r,t) 

To solve the general case we take the following two assumptions: 
[i) v = 0, A = A(r, t); (ii)v = v(r, t), A = 0. 
(i) When v = 0, A = X(r,t) 

EFE (6) implies that A = In fR' , where f(r) is an arbitrary function of 
r and R is an arbitrary function of coordinates r and t. Also R' ^ 0, which 
implies that R ^ R(t). 

Here arises two cases i.e. either R = R(r) or R = R(r, t). For R = R(r), 
we can have p and p by replacing A in Eqs. (3) and (4) respectively. 

9 k R? f 2 R 2 RR'f 3 



But Eq.(5) gives 



P ~ k { R 2 + pR?'' 



f 

P kRR'P' 
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By comparing the two values of p, we obtain 



where / is an integration constant. The resulting metric becomes 

This turns out a class of spherically symmetric static spacetimes. For / = 
-y== and I — 1, it reduces to Einstein metric. From the above equations we 

obtain p — V = ~ l ~ which implies that p + 3p = 0. 

(ii)When v = u(r, t), A = 0, EFE (6) will reduce to v — In g R , where 
g(t) is an arbitrary function of t and R is an arbitrary function of coordinates 
r and t. Also 0, which implies that R ^ R(r). This gives that either 
R = R(t),R = R(r,t). For R = R(t), we have the spacetime similar to 
Eq.(16). The quantities p and p also turn out to be the same. 

For R = R(r,t), the solutions need to be investigated. 

4 Kinematic of the Velocity Field 

The spherically symmetric solutions can be classified according to their kine- 
matical properties [6]. The rotation is given by 

Uab = U [a ;b]+ U[aU b ]. (30) 

Tha acceleration can be found by 

U a = U a - b U b . (31) 

For the expansion we have 

= u%. (32) 
The components of the shear-tensor are given by 

Oab = «(a;6)+ «(a U b ) - ^Qh ab , (33) 
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where h ab = g ab — u a u b is the projection operator. The square brackets 

denote antisymmetrization and the round brackets indicate symmetrization. 
The shear invariant is given as a a b<7 ab . 

Now we find all the above quantities for the solutions obtained. The 

rotation and the acceleration are zero for all the solutions. The expansion, 
for the first solution, is 

e = (^)*[££^ + -] (34) 

v t ' l 2t(c-t) V y ' 

The components of the shear-tensor are given by 



2 c - t) i 1 put 3 - c 



3 V t ' l t 2t(c-ty 



c 2 {l-(-— 1 (-) 2 + (^) 2 -} 

6 C Z ri 
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(35) 



1 9 .c — t) . i r p«t 3 — c 1, . . 

, 22= _ (2( _i) 3[ |___], (36) 

o- 33 = sin 2 #cr 22 . (37) 
For the second solution, the expansion factor is 

The components of the shear-tensor are given by 

(l-a 2 t 2 )5 / x 

a 22 = g— A (40) 

o- 33 = sin 2 #cr 22 . (41) 



For the third solution, we have 



The components of the shear-tensor are given by 



e = + ( 42 ) 



°" = 5<7 " £ M ' (43) 



•* = 5**<s-7Va < 44 > 

a 3 3 = sin 2 9a 2 2- (45) 

The shear invariant turns out to be 3 for all the solutions. 

The rate of change of expansion with respect to proper time is given by 
Raychaudhuri's equation [12] 

^ = -^e 2 - a ab a ah + u ab u a u b - R ab u a u» (46) 
dr 3 

We evaluate it only for the second solution as it is simple to understand. For 
this solution it becomes 

d& 2 ,1 + aH 2 . 

Tr = ^^^ ) (47) 



5 SUMMARY 

This section presents a brief summary of the results obtained in the previous 
section. We shall discuss these results and will comment on possible future 
directions of development. 

We have been able to solve partially the EFEs for the classes of non- 
static spherically symmetric spacetimes using different equations of state for 
perfect fluid. These have been classified into three categories. The summary 
of each case is given below: 

1. R = R(t), A = A(t) 

In this case, the dust solution gives 



ds 2 



t 



-dt 2 - 



c-t 

The stiff matter solution becomes 

a 2 t 2 



c 2 {i-(-V)^^ _1 (^ + (V^7> n2 

L C V ri 



ds z 



1 - a 2 t 2 



dt z - dr l - t 2 dn 2 



dr 2 - t 2 dtt 2 
(48) 

(49) 



2. R = R(t), A = A(r,t) 
The stiff matter solution gives 
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ds 2 = [y{t)]- l dt 2 - z 2 (r, t)dr 2 - t 2 dtl 2 , (50) 

where y and z are already given. 
3. R = R(r,t) 

In this case we obtain two solutions. The first solution, infact, turns out 
a class of spherically symmetric static spacetimes. The other (stiff matter) 
solution is exactly the same as Eq.(49). 

The non-static spherically symmetric solutions with equations of state 
split up into three classes of solutions. In the first case we obtain two different 
solutions. One of these two becomes the dust solution and the other becomes 
a stiff matter solution. The pressure and the energy density are positive every 
where. In case two we have a stiff matter solution only. However, it is much 
difficult to comprehend it. In the 3rd case we obtain two solutions. The first 
solution becomes a class of spherically symmetric static spacetimes depending 
upon the arbitrary function f. If we choose / = ^^_ r2 , it reduces to Einstein 
spacetime. The energy density is positive everywhere while the pressure is 
negative. The other (stiff matter) coincides with one of the solutions in the 
first case. It is interesting to note that all the non-static solutions which have 
been found contain shear which are a very few solutions in the literature. 

Finally, we discuss the behaviour of the rate of change of expansion for 
one solution given by Eq.(49). We see from Eq.(47) that the rate will never be 
positive but this will always be negative. As the time t tends to zero, the rate 
approaches to — oo and as t goes to oo, the rate tends to zero. This shows 
that the spacetime is contracting or collapsing and the flux gets focussed 
along the proper time. Similarly, the other solutions can be discussed. 

We have tried to obtain non-static spherically symmetric solutions for 
some particular classes and partial solutions have been attempted in these 
three classes. Thus a total of three solutions have been obtained. To obtain 
new solutions one has to solve the remaining cases of these classes. Then 
we can attempt the general solution of non-static spherically symmetric so- 
lutions. 
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